Structural effects on the spin-state transition in epitaxially strained LaCoOa films 
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Using density functional theory within the LSDA + U method, we investigate the effect of strain 
on the spin state and magnetic ordering in perovskite lanthanum cobaltite, LaCoOs. We show 
that, while strain-induced changes in lattice parameters are insufficient to stabilize a non-zero spin 
state, additional heteroepitaxial symmetry constraints - in particular the suppression of octahedral 
rotations - stabilize a ferromagnetic intermediate-spin state. By comparing with experimental data 
for the bulk material, we calculate an upper bound on the Hubbard U value, and describe the role 
that the on-site Coulomb interaction plays in determining the spin-state configuration. 
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I. INTRODUCTION 

The desire to control magnetism with external per- 
turbations other than magnetic fields has motivated 
much recent research on the s train - and electric field- 
response of magnetic materials) ^ ! ^ * '^ ! Such control might 
enable future Mottronic applications, in which small 
external perturbations drive transitions between com- 
peting electronic, orbital, charge and spin orderings, 
causing drastic changes in properties. Thin film het- 
erostructures containing transition metal oxides are prov- 
ing particularly promising in this emerging field; for 
example, electric-field switching of magnetization has 
been achieved through exchange-bias coupling of ferro- 
magnetic Coo.gFeo.i to multiferroic BiFeOs, and sub- 
strate induced strain has been used to tune magnetic 
interaction^ in magneto elastic c omposites (see Ref.|B]for 
a review). Recent report J^MESl of a substrate-dependent 
spin state in epitaxial films of LaCoOa are of partic- 
ular interest since they suggest a route to switching 
magnetism off (low spin diamagnetic (f' Co^+) and on 
(intermediate- or high-spin Co'^+). 

Lanthanum cobaltite is a rhombohcdral (i?3c) per- 
ovskite that has been of continued interest for the last 
half-century, due in part to the many magnetic phase 
transitions that occur as a function of t emperature, pres- 
sure and chemical doping .^^E^E^ISIil! These transitions 
are a consequence of the competing crystal-field split- 
ting energy (Aqf), Hund's exchange energy (Jh) and 
d-orbital valence bandwidth {W), which are similar in 
magnitude, resulting in low-, intermediate- or high-spin 
Co'^"'", depending on the details of the system. In the 
ground state (T = K), LaCoOa is a diamagnetic insu- 
lator with a low-spin (S'=0, ^Ig^g) Co'^"'" configuration. 
It is thermally excited to a paramagnetic intermediate- 
{S=l, t2geg) or high-spin (5=2, tig^l) semiconducting 
state above approximately 95 K.^^ The nature of this 
spin-state transition is still under debate: inelastic neu- 
tron scattering, x-ray absorption spectroscopy (XAS) 
and magnetic circular dichroism experiments'^ suggest a 
first-order transition to the high-spin (HS) state, while 
other x-ray photoemission (XPS) and XAS spectra in 
addition to electron energy loss (EELS) spectroscopy 



suggest the intermediate-spin (IS) state.'^^EHUnill] gim- 
ilarly, Hartree-Fock clustei'221 and full-potential DFT 
calculations^^ suggest the HS state is more stable than 
the IS state, whil e oth er LSDA -t- U calculations obtain 
the reverse result .'^^HH 

In contrast to the bulk behavior, recent studies on 
LaCoOa thin films report ferromagnetism with a field- 
cooled magnetization of 0.37 /is /Co ion on a sub- 
strate that causes 1.84% tensile strain.*^ Although 
ferromagnetic hysteresis loops have been recorded by 
several groups,l^ ^ * ^ * ^° * ^^ l it remains unclear experimen- 
tally whether such magnetism is an intrinsic feature of 
strained LaCoOs, or whether it arises from sample off- 
stoichiometry (ferromagnetism induced by hole doping 
is observed in bulk Sr-rich LaCoOs samples) or is 
a surface effect resulting from the change in coordina- 
tion of surface Co ions (also recently reported in bulk 
samples) ^ 

In this work, we use ab-initio calculations based on 
density functional theory (DFT) to show that epitaxial 
strain can indeed drive a spin-state transition to a fer- 
romagnetic state in stoichiometeric LaCoOa. The tran- 
sition is not caused, however, by strain-induced changes 
of the lattice constants, but rather relies on interface- 
induced changes in the tilt pattern of the CoOg octahe- 
dra. 



II. COMPUTATIONAL DETAILS 

We use the the projector augmented plane wave (PAW) 
method of DFTj^SSl as implemented in the V ienna Ab Ini- 
tio Simulation Package (vASP) code .1221221 To accurately 
describe the exchange and correlation, we use the spher- 
ically averaged form of the rotationally invariant local 
spin d ensity approximation + Hubbard U (LSDA+C/) 
methocPUaai^ith 

one effective Hubbard parameter UeS = 
U — J, and treat the double counting term within the 
fully localized hmit. We use the supphed VASP PAW 
pseudopotentials (La_S, Co, 0_s) with the 5p'^5d^6s^ 
valence configuration for La, 4s^3d* for Co, and 2s^2p'^ 
for O. Other technical details include a plane wave energy 
cutoff of 550 eV, a 7 x 7 x 7 fc-point grid to sample the 
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Brillouin zone, and the tetrahedron method with Blochl 
corrections^ and an 11x11x11 fc-point grid to calculate 
the densities of states. 



III. BULK LaCoOs 
A. Correlation Effects 

The LSDA+ Hubbard U approach has been success- 
ful in treating static correlations in transition metal ox- 
ides; however the selection of an appropriate UeS is rarely 
straightforward, and a number of methods exist for de- 
termining suitable values. These include experimental 
measurement from p hotoemi ssion spectroscopy,'^^! self- 
consistent calculation d^^l'^^ l ^ and educated guesswork. 
LaCoOa represents a particularly difficult case because 
of the strong dependence of orbital occupation - which 
affects the polarizability and screening, and in turn the 
UeS - on pressure and strain. 

Indeed, XPS experiments have found the d-d Coulomb 
correlations to range from 3.5 to 7.5 eVi-^^^ depending 
on the structural details of the samples. Previous single- 
site and two-site configuration-interaction cluster calcu- 
lations have obtained values from 4 to 5.5 eV."'^^ On the 
other hand, recent first-principles calculations suggest 
that the spin state is independent of the choice and 
values as large as 9 eV have been used to study temper- 
ature dependent spin-state transitions.'^^! Each method 
does agree however that the d-d electron repulsion for 
the low- and intermediate-spin states is approximately 
the same. 

Due to these discrepancies in the literature, in the first 
part of this study we revisit the effects of electron re- 
pulsion on the spin-state and orbital occupation in bulk 
LaCoOa. We begin by calculating the critical U value 
that induces a spin state transition in the bulk rhombo- 
hedral structure by calculating the total energies for each 
spin state as a function of Ucff. Since the bulk zero kelvin 
state is known to be low spin, this requirement then pro- 
vides a bound on allowable Ucs values for studying the 
experimentally observed ferromagnetic thin films. 

Here we use the low temperature experimental struc- 
tural parameters found in Ref. [Tni (a = 5.275 A, a = 
61.01°) and fully relax the internal coordinates with 
Ucff = eV until the forces are less than 1 meV A~^. 
Within the LSDA, we find the correct ground state struc- 
ture: a diamagnetic insulator with a 0.45 eV band gap, 
which is close to the measured optical gap.*^ Our cal- 
culated energies and magnetic moments as a function of 
UcB are shown in Figure [l] The most striking finding is 
that the experimentally observed S = ground state is 
only stable for Ucs values less than 4.0 eV. Therefore we 
regard 4.0 eV as an upper bound on UeS for LaCoOa. 

In Figure [ijjinset) the relative energies of the diamag- 
netic S" = 0, IS = 1, and HS 5 = 2 states are also 
shown. [For these comparison calculations we impose fer- 
romagnetic (FM) order in the IS and HS Co sublattice 
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FIG. 1: (Color online) Calculated magnetic moment as a 
function of Uc« in bulk rhombohedral LaCoOa with the cor- 
responding orbital energy diagrams, (inset) Relative energies 
per formula unit of the 5 = 0, S=l and S=2 spin states. 



so that we can fix the total spin moment.] As expected, 
as Ucff is increased, spin pairing in the manifold be- 
comes less favorable as the energy gain from the Hund's 
exchange energy exceeds the energy cost in creating a 
singlet state and thereby reduces the relative energies of 
higher spin states. Such correlation-induced spin-state 
transitions, in which higher Ueff values induce states with 
higher magnetic moments, have previously been found in 
a number of other transition metal compounds. ""^ 4i42] 
Co'^"'" is a particularly interesting case, however, because 
the low-spin state is non-magnetic and so the spin-state 
transition increases the magnetic moment to a finite value 
from an initial value of zero. Interestingly, at the transi- 
tion to the 5=1 state, we find Ucfi /W—Q.27 , which is 
low compared to most moderately- or strongly-correlated 
magnets; in addition, W is largely independent of Ucft 
(not shown). For all Ucff values we find that the HS 
state is more than 1.0 eV higher in energy than the IS or 
LS states. 

The critical UaS value we have determined with this ap- 
proach is in good agreement with x-ray photoemission ex- 
periments reported in Ref. and recent first-principles 
calculations in Ref. |43l Therefore, for the remainder of 
this study, we strictly use values of [/cff < 4.0 eV unless 
noted otherwise. The strong dependence of the ground- 
state spin configuration on [/off partly explains the incon- 
sistencies between different first-principles calculations in 
describing the evolution of the spin-state transition. 



B. Electronic Structure 

Before investigating the effects of strain on the mag- 
netic behavior we describe the nature of the unusual in- 
termediate S = 1 state of LaCoOa compared to the dia- 
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FIG. 2: (Color online) Spin- and orbital-resolved densities of 
states for S = (upper) and 5=1 (lower) rhombohedral 
LaCoOa {UcB ~ 3.0 eV). In the low-spin state, a diamagnetic 
insulator is found, while a half-metallic ground state is found 
with a local magnetic moment of 1.8 hb per Co atom in the 
intermediate spin state. 



magnetic S = state. In the molecular cluster limit [Fig- 
ure [ij inset)], when AcF > Jr the low-spin configuration 
is favored, while when A^p < Jjj the high-spin state 
dominates due to the gain in exchange energy from the 
parallel alignment of spins; the intermediate spin state 
might be expected when Acf ~ Jh- Furthermore, as evi- 
dent in Figure[2] hybridization and covalency between the 
O 2p and Co 3d states causes dramatic deviations from 
the simple molecular cluster picture by causing strong 
broadening of the bands; in particular the Co Cg orbitals 
span more than 11 eV in energy. 

In Figure [2] we show our calculated electronic densities 
of states for the S = and constrained S = 1 ferromag- 
netically ordered LaCoOa at the experimental lattice pa- 
rameter with a UeB = 3.0 eV. In the low-spin state, the 
Co t2g manifold is fully occupied with triply degenerate 
dxy, dxz and dyz orbitals; the valence band is formed by a 
mixture of these Co states and O 2p orbitals. The doubly 
degenerate d^^-y^ and dz^ orbitals which make up the Cg 



manifold form the conduction band. 

In the intermediate-spin state, on the other hand, 
broadening of the band widths causes the majority up- 
spin eg band to extend below the Fermi energy, and the 
minority t2g band to extend above. The hole in the 
<2g manifold consists of a superposition of the minority 
-^{d^y + dyz + dxz) orbitals, and the "missing" electron 
occupying the lower part of the majority Cg band is in 
a ;^(c^z2 -f dx'i-y^) state. This behavior causes our cal- 
culated ferromagnetic intermediate-spin state to be half- 
metallic. Note, however, that bulk LaCoOs is in fact an 
insulator up to room temperatur^^ we will return to 
this discrepancy later. 



IV. STRAINED LaCoOs 

Spin-state transitions in bulk LaCoOa are known to 
occur as a function of unit cell volume,'^ where the 
tendency for an excited spin state is favored for larger 
unit cell volumes due to competition between the en- 
ergy penalty in forming a singlet state and the gain 
in Hund's exchange energy favoring ferromagnetic align- 
ment of spins. This is evidenced by the fact that high 
spin Co'^^ has a larger ionic r adius (0.61 A), compared 
to the low spin Co3+ (0.55 A) .Esl Indeed, this IS consis- 
tent with our calculations for the ideal cubic pervoskite 
where the LSDA-|-[/ equilibrium volume of the interme- 
diate state is approximately 2% larger than that of the 
low spin configuration. 

In this section we examine likely crystallographic struc- 
tural distortions in thin films to determine whether they 
cause a spin-state crossover. In general, heteroepitax- 
ial strain from coherent growth on a substrate with a 
mismatched lattice constant can modify the structure by 
changing the lattice parameters, symmetry or chemistry 
at the interface. Therefore we explore whether compu- 
tations that incorporate changes in Co-O-Co bond an- 
gles, Co-0 bond lengths, unit cell volume, or combina- 
tions of these effects, are able to reproduce the ferromag- 
netic state which is observed experimentally in thin film 
LaCoOs. 



A. Effect of Changes in Lattice Parameter 

The low temperature rhombohedral structure of 
LaCoOa belongs to the (a^a^a^) Glazer tilt system, 
in which successive octahedra rotate in opposite senses 
along each crystallographic direction. The Co-O-Co 
bond angle is approximately 166°. The importance of 
such octahedral rotations on the electronic properties of 
thin film perovskite oxides has been the subject of many 
recent reports, pa rticul arly in the context of their effect 
on ferroelectricity .SSllH effg^^ on magnetic properties 
is also Hkely, because changes in Co-O-Co bond angles 
can strongly affect the magnetic superexchange inter- 
actions. In this section we investigate whether strain- 
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FIG. 3: (Color online) (a) Energy difference between the 
diamagnetic low-spin and ferromagnetic intermediate-spin 
states, Ai5, as a function of epitaxial strain applied in the 
pseudo-cubic (ps) (100) plane of the rhombohedral structure 
(b) shown relative to the cubic perovskite. (c) The relaxed 
mean Co-0 bond lengths, bond differences &d between the 
long- and short bonds in the CoOe octahedra, and the O- 
Co-O bond angles for C/efi = 3.5eV. The lines are a guide for 
the eye. 



induced changes in these bond angles are sufficient to 
stabilize the intermediate-spin state. 

In Figure [sjj a), we plot our calculated energy differ- 
ence between low- and intermediate-spin states with fer- 
romagnetic order as a function of strain applied to the 
pseudo-cubic (100) plane with respect to the LSDA equi- 
librium volume. Values are shown for two J7off values, 3.5 
and 4.5 eV, chosen to be above and below the critical U^s 
value of 4.0 eV that we established in the previous sec- 
tion. (The high spin state is not shown, since it is ~1 eV 
higher in energy at all strain values.) At each in-plane 
strain value we adjust the out-of-plane lattice parameter 
and rhombohedral angle to maintain the bulk equilibrium 
volume, then fully relax the atomic positions. We find 
that, for i7eff below our calculated critical value, the low- 
spin ground state is stable up to strains of 4%; therefore 
we predict that strain-induced changes in lattice parame- 
ters alone are insufficient to cause a spin-state transition 
in LaCoOs up to reasonable strain values. Even our un- 



physically large ?7cff of 4.5 eV does not induce a transi- 
tion to the intermediate spin state until just beyond 4% 
compressive strain. The half-metallic state remains sta- 
ble until strain values above 3% (J7cff — 3.5 eV) and 1% 
(?7eff = 4.5 eV), when the electronic structure becomes 
fully metallic. 

To understand the absence of spin-state transition with 
strain, we plot in Figure [3|^c) the evolution of the mean 
Co-0 bond length and the mean Co-O-Co bond angle 
as a function of strain for IJaS = 3.5 eV. With either com- 
pressive or tensile strain, the average Co-0 bond length 
increases from the equilibrium (zero-strain) value. All 
bond lengths increase uniformly, however, such that the 
CoOe octahedra remains perfectly octahedral, and the 
ideal Oh crystal field is maintained. This is supported by 
the bond length differences between the long- and short- 
Co-0 bonds ((5c?) in the CoOg octahedra. We find signif- 
icant changes in the Co-O-Co bond angles, particularly 
for tensile strain. In many magnetic perovskites, such 
large changes in b ond a ngles are sufficient to change the 
magnetic ordering .I^^III! believe that the absence of 
spin crossover in this case is due to the exceptionally 
broad bandwidth of the Co orbitals, which reduces 
the exchange energy gain from spin polarization. 



B. Effect of Octahedral Rotations and Distortions 

Next we isolate the influence of these octahedral rota- 
tions by manually disabling them while applying strain. 
Our motivation is two-fold. First, there is experimental 
evidence suggesting that LaCoOs gro ws in s uch a 'cube- 
on-cube' manner on many substrates! ''1^1^° ! Second, dis- 
abling the octahedral rotations causes the system to re- 
spond to strain by changing the local symmetry around 
the Co ion; therefore we can examine the influence of the 
crystal field on the spin-state transition. 

The no-rotations constraint is imposed by using a five 
atom unit cell which prohibits rotations by symmetry; 
as a side-effect this also imposes ferromagnetic ordering. 
(We later examine if this is indeed the preferred magnetic 
ordering.) We begin by setting the in-plane pseudo-cubic 
lattice parameter (a) to that of the experimental sub- 
strate (LaA103)o.3(Sr2AlTa06)o.7 (LSAT) value (3.87 A) 
with the optimized LSDA-|-?7 volume (c/a — 0.955), and 
relax the internal coordinates and out-of-plane (c) lat- 
tice constant; the resulting structure is diamagnetic and 
270 meV higher in energy than the LSDA equilibrium 
-R3c LS structure. 

Next, we apply uniaxial strain by varying the c/a ra- 
tio with a fixed to the experimental LSAT value, and 
show the resulting calculated magnetic moment in Fig- 
ure |4] Our main finding is that, at the experimental c/a 
ratio, the IS state is lower in energy than the LS, and 
therefore we predict that the Co ions should be mag- 
netic. The origin of the stabilization of the intermediate 
state is the lifting of the octahedral crystal field by the 
tetragonal symmetry adopted when the octahedral rota- 
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tetragonal distortion (c/a) 

FIG. 4: (Color online) Calculated magnetic moment as a func- 
tion of c/a for tetragonal LaCoOa with a fixed to the exper- 
imental LSAT lattice parameter (3.87 A). The dashed line 
(filled circle) indicates the experimental (LSDA equilibrium) 
c/a ratio, (inset) Schematic representation of Co displace- 
ment within the equatorial octahedral plane (arrow). 
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TABLE I: Total energy differences (in meV) for various mag- 
netic orderings within the tetragonal crystal structure relative 
to the five atom ferromagnetic unit cell. Values are given at 
the experimental and optimized LSDA+(7 c/a ratios. Cal- 
culated magnetic moments per Co atom in are given in 
parentheses. 



tions are disabled. Small (< 1%) uniaxial expansion of 
the out-of-plane lattice constant modifies the tetragonal 
crystal field splitting sufficiently to favor occupation of 
the Cg manifold. In contrast, when octahedral rotations 
are allowed, strain is accommodated through changes in 
the rotation angles rather than through modification of 
the local bond lengths around the Co ions; the local crys- 
tal field splitting is therefore largely unchanged and the 
diamagnetic state remains stable. 

Next we investigate whether the ferromagnetic order 
imposed so far for computational convenience is indeed 
the lowest energy magnetic ordering by comparing its 
energy with those of the A- type antiferromagnetic (a- 
AFm), and G-type antiferromagnetic (g-AFM) orderings 
for this structure. The total energies for each structure 
are shown in Table [l] for both the experimental and opti- 
mized LSDA+C/ c/a values. The energies are given rela- 
tive to the FM single unit cell configuration. We find that 
the G-AFM and A- AFM structures are 180 and 375 meV 
per formula unit higher in energy than the ferromagneti- 
cally ordered intermediate spin state at the experimental 
c/a ratio. This is consistent with the experimentally ob- 
served ferromagnetism. 

Intriguingly, the structure with disabled octahedral ro- 
tations relaxes to a polar space group (C2/c) with the 
Co^"*" ion moving 0.06 A off-center in the ah plane. By 
summing the formal ionic charges multiplied by their dis- 



placements from their centrosymmetric positions we find 
an in-plane polarization of 17.5 /LtC/cm^ at the experi- 
mental c/a = 0.982 ratio. Note that since our overall 
electronic ground state is metallic, we are unable to eval- 
uate the electronic contribution to the polarization using 
the standard Berry's phase approach. Indeed the onset 
of ferroelectric polarization caused by the disabling of oc- 
tahedral rotations in perovskite oxides has been noted in 
a number of calculations, and is believed to result from 
the o ff-cen tering of ions to maintain a favorable bond 
order.EsEZlUowever, since our electronic structure is over- 
all metallic, we cannot predict ferroelectric behavior. 



V. DISCUSSION 

While our finding of ferromagnetism is consistent with 
recent experimental reports, there are some important 
differences between our computations and the experi- 
mental observations. First, and analogous to the bulk 
intermediate-spin case, our calculated tetragonal struc- 
ture is half-metallic, with a broad majority spin O 2p - 
Co Cg band crossing the Fermi level; experimentally the 
ferromagnetic films are found to be insulating. In addi- 
tion, the sizes of most measured magnetic moments are 
an order of magnitude smaller than our calculated value. 
Recent x-ray magnetic circular dichroism (XMCD) ex- 
periments on thin films, however, find a local Co mo- 
ment of 1.2 Up, which is in better agreement with our 
calculations 

In this last section we attempt to reconcile our finding 
of a strain-induced half-metallic ferromagnetic arrange- 
ment with experimental reports of insulating LaCoOa on 
LSAT. In particular, we explore likely Jahn- Teller dis- 
tortions and orbital orderings which are known to allow 
both ferromagnetism and insulating behavior in oxides.-^ 
We also examine the effect of including spin-orbit interac- 
tions in our calculations, since these couplings can make 
significant contributions to determining the orbital occu- 
pation in many transition metal oxides. 

The possibility of an orbitally ordered state in LaCoOa 
was suggested previously in Ref. 1231 and unrestricted 
Hartree-Fock calculations on similar material^^ found 
that small Jahn- Teller structural distortions can stabilize 
an insulating state. Cooperative Jahn- Teller distortions 
(ranging from 1 to 6% from low to room temperature) 
have indeed been demonstrated in LaCoOa with various 
techniques including high-resolution x-ray diffraction,'^ 
Raman scattering,^** and neutron diffraction.'^SM] j^gf. 
1531 obtained a monoclinic structure with 12 /a symmetry, 
consistent with the A- AFM ordering seen in LaMnOa. 

Although the type-A monoclinic structure was investi- 
gated previously in Ref. dSJ we revisit the possible orbital 
ordering available in strained LaCoOs by imposing from 
1-6% Jahn- Teller structural distortions of the CoOg- In 
the same way, we study type-d ordering (space group 
PA/mbm), which has uniform orbital occupation along 
the c-direction and alternating orbital occupation in the 
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FIG. 5: The type-d Jahn- Teller distorted structure is shown 
with the possible eg orbital ordering configuration. The or- 
bital degeneracy is also split at each of the two sites, according 
to the orientation of the elongation of the oxygen octahedra 
with respect to the c lattice parameter. 

a&-plane, consistent with overall ferromagnetic superex- 
change (Figure [5| . This 2D antiferrodistortive behavior 
allows for (1^2 _y2 and orbital ordering due to the spht- 
ting of the eg degeneracy from the structural distortions 
of the oxygen octahedra. Similarly the threefold degen- 
erate t2g orbitals also spht into a two-fold degeneracy, 
that is lower in energy than the single non-degenerate 
state. In all cases examined no Cg orbital order is ob- 
served, in spite of the imposed octahedral structural dis- 
tortions. Due to the degeneracy of the Co Cg states at 
the Fermi level, although the density of states is reduced, 
a half-metallic ferromagnetic state persists. In contrast, 
the higher energy G-AFM structure considered earlier is 
insulating. 

Finally, we examine the effect of spin-orbit interactions 
in LaCoOs, since such coupling of the magnetic spin de- 
grees of freedom to the lattice can change the relative 
level splitting and degeneracy of spin a nd or bital ground 
states in the Sd transition metal oxides .'^^'^ For the dia- 
magnetic S — Q configuration we find that spin-orbit in- 
teraction has a negligible effect on the electronic struc- 
ture as expected for a filled t2g manifold. On the other 
hand, in a recent theoretical study^'^ it was suggested 
the S = 1 state does show significant changes in the elec- 
tronic structure near the Fermi level when spin orbit in- 
teractions are included. In contrast to those results which 
used a generalized-gradient approximation (GGA-f f7) for 
the exchange-correlation functional, we do not find sig- 
nificant deviations in the electronic structure when spin- 
orbit coupling is included in the calculations. This result 
is likely due to the difference in the description of our 
starting intermediate spin-state configurations: we begin 
with a half-metallic ground state, whereas in Ref. |35J a 
fully metallic state is found with a high density of states 
at the Fermi level. When spin-orbit coupling is included 
in the calculation, the large peak feature is naturally split 
into a doublet. 

The origin of the two inconsistencies - our half- metallic 



rather than insulating ground state, and larger mag- 
netic moment per Co ion compared with experiment - 
might lie in the difficulties associated with producing 
and characterizing high quality, uniform transition metal 
oxide films, or from a failure of the LSDA+C/ method 
to fully describe the complex orbital physics. Future 
theoretical investigations should consider more sophis- 
ticated methods such as dynamical mean field theory, 
in which dynamical correlations (spin ffuctuations) can 
in principle be treated explicitly. Indeed, recent in- 
elastic neutron scattering measurements report a small 
Jahn- Tell er dis tortion which has short-range dynamical 
character l^^l^^ l This dynamic Jahn- Teller effect is consis- 
tent with a proposed vibronic e^-0-e° superexchange be- 
tween intermediate-spin Co atoms and would allow for 
fluctuations of AFM exchange which should reduce the 
magnetic moment. On the experimental front, our calcu- 
lations suggest that more detailed characterization of the 
local electronic and structural properties will be invalu- 
able in understanding and exploiting the spin behavior 
of LaCoOa films. 



VI. CONCLUSIONS 

In summary, by comparing our calculated LSDA+C/ 
spin state of bulk LaCoOa with the measured low tem- 
perature behavior, we have determined a critical upper 
bound of 4 eV on the Hubbard U parameter for LSDA+C/ 
calculations for LaCoOa. Using our critical U value, 
we have established that strain-induced changes in lat- 
tice parameters are insufficient to cause transitions to 
finite magnetic moment spin states at reasonable val- 
ues of strain. Instead, if the cooperative octahedral tilt- 
ings and rotations are deactivated, intermediate-spin lo- 
cal moments are stabilized on the Co ions at small strain 
values, and these order ferromagnetically. Our results 
suggest a possible route to dynamically controlling mag- 
netism using an electric field, in superlattices of LaCoOa 
with a piezoelectric material. 
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